Electrical transport properties of stable single-atom contacts of zinc 
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We report low-temperature measurements of the electrical transport properties of few-atom con- 
tacts of the superconducting metal zinc, arranged with lithographically fabricated mechanically 
controllable breakj unctions (MCB). The conductance histogram shows several narrow peaks, not 
regularly observed for multivalent metals. The first peak - corresponding to the single-atom contact 
- is located slightly below one conductance quantum and is split into two subpeaks, possibly indi- 
cating two preferred configurations of the single-atom contact. The current-voltage characteristics 
in the superconducting state show nonlinearities due to multiple Andreev reflections (MAR). In 
addition we find nonlinearities in the current- voltage characteristics that we attribute to electronic- 
interaction effects due to the granular structure of the film. Some of these latter findings are in 
qualitative agreement with a prediction by Avishai et ali 
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For revealing and understanding the electronic trans- 
port properties of atomic scale circuits, single- or few- 
atom contacts are used as model systems. In this 
limit the electronic conduction can be regarded as a 
fully quantum-mechanical scattering and wave-matching 
problem, i.e. the conductance can be attributed to inde- 
pendent electronic modes, nick-named "conduction chan- 
nels" . An atomic-size contact between two metallic elec- 
trodes can accommodate only a small number of conduc- 
tion channels. The transport properties are thus fully de- 
scribed by a set {t„} = {n, T2, ...tk} of transmission co- 
efficients which depends both on the chemical properties 
of the atoms forming the contact and on their geometrical 
arrangement. Experimentally, contacts consisting of even 
a single atom have been obtained using both scanning- 
tunnel microscopy and breakjunction techniques 2 . The 

N 

total transmission D— t; of a contact is deduced from 

i=l 

its conductance G measured in the normal state, using 
the Landauer formula G — GqD where Go = 2e 2 /h is 
the conductance quantum 3 ". 

Experiments on a large ensemble of metallic contacts 
have demonstrated the statistical tendency of atomic-size 
contacts to adopt element-specific preferred values of con- 
ductance. The actual preferred values depend not only on 
the metal under investigation but also on the experimen- 
tal conditions. However, for many metals, and in particu- 
lar 'simple' ones (like Na, Au...) which in bulk form good 
'free electron' metals, the smallest contacts have a con- 
ductance G close to Go (Ref. 2 ). Statistical examinations 
of Al point contacts at low temperatures yield preferred 
values of conductance at G = 0.8 Go, 1.9 Go, 3.2 Go and 
4.5 Go (Refi^), indicating that single-atom contacts of 
Al have a typical conductance slightly below the conduc- 
tance quantum. 

However, it has been shown by the analysis of current- 
voltage (IV) characteristics in the superconducting state 
that Al single-atom contacts accommodate in general 



three transport channels, the transmissions of which 
add up to a value around 1 (RefpS.). These find- 
ings are in agreement with measurements of the shot 
noise&i, conductance fluctuations^, thermopower— and 
the supercurrentiifl in such contacts. A quantum chemical 
modeii^ which links the r n to the chemical valence and 
the atomic arrangement of the region around the central 
atom of the constriction, has been put forward. This 
model predicts that monovalent metals would transmit a 
single channel while sp-like metals as e.g. Al or Pb up 
to four channels due to the contributions of the (unoccu- 
pied) p orbitals. So far Al appears to be the only metal 
with the property to show peaks in the histogram close 
to multiples of Go although this does not correspond to 
quantized transmission values. 

According to simplistic arguments and ignoring the 
band structure, divalent elements such as the alkaline- 
earth elements or the IIB subgroup elements Cd, Zn and 
Hg, should be insulating, since they have a completely 
filled outer s-shcll. This apparent discrepancy to the 
experimental findings indicates that additional orbitals 
besides the s-orbitals contribute to the electronic con- 
ductance. Tight-binding calculations reveal that in the 
case of Zn the 4p orbitals are the most relevant ones for 
a correct description of the bulk band structure (Refii). 
With these considerations, single-atom contacts of Zn are 
expected to have very similar transport properties as Al 
single-atom contacts, i.e. up to four channels with a to- 
tal conductance of about 1 Go^. The distribution of the 
transmission coefficients, however, may deviate, since it 
depends on the exact atomic configuration. Since Zn 
crystallizes in a hexagonal, slightly distorted hep struc- 
ture while Al is a fee crystal, the actual atomic arrange- 
ment of the point contacts might be different. Moreover, 
due to the "loss" of many neighbors, a single Zn atom in 
a contact may behave more atomic-like. 

Here, we report the fabrication and low-temperature 
electronic transport properties of lithographic MCBs of 
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FIG. 1: Electron micrograph of sample No. 1 under an 
inclination angle of 70°, taken after measurement. The poly- 
crystalline structure of the Zn film is clearly visible. A re- 
maining spacing of the bridge arms of several nm and a fine 
nanowire between them is observed. 



Zn. This work represents the first study of atomic-size 
contacts of Zn. The measurements in the normal state 
demonstrates the tendency of the atomic-size contacts 
to adopt preferred conductance values. The analysis of 
the current-voltage characteristics reveals superconduct- 
ing properties of the nanowires as well as indications 
of electron-electron interactions. We show data of four, 
nominally identical samples, labeled No. 1 to 4. 

The samples are fabricated along the lines of Ref^. 
Our samples are 2 /im long, 200 nm thick suspended Zn 
nanobridges, with a 200 nm x 100 nm constriction in the 
center (cf. Fig.^l. The thickness of the bronze substrates 
is w 0.3 mm, the distance between the counter-supports 
of the bending mechanism is 10 mm. With these sample 
parameters and the geometry of the mechanics we achieve 
reduction ratios (between the motion of the pushing rod 
and the displacement of the bridge arms) of the order of 
5000:1 (see below). 

Special care has to be given to the evaporation of Zn. 
Because of its low boiling temperature and high vapor 
pressure, Zn desorbs easily from the substrate if no crys- 
tallization seeds are provided. If Zn is evaporated at 
room temperature, it grows in islands of micron-size with 
the hexagonal axis preferably perpendicular to the sub- 
strate plane. Since this structure is not suitable for the 
formation of atomic-size contacts, we take special care 
to reduce the grain size. We first deposit a 1 nm thick 
seed layer of Ge onto the substrate which is cooled by 
liquid nitrogen to a temperature of w — 50°C. At lower 
substrate temperatures the different thermal expansions 
of the metallic substrate, the Zn itself, the polyimide sac- 
rificial layer and the polymer lithographic mask provoke 
cracks in either the mask or the Zn layer on warming 
up to room temperature. The metal is then evaporated 
from a home-made boron-nitride crucible heated by a 
tungsten foil at a pressure of m 10~ 6 mbar and a rate of 
w 0.2 nm/s at the same substrate temperature. The pa- 



rameters are chosen empirically such that the grain size 
of the Zn layer does not exceed « 200 nm and the ad- 
hesion is good enough for the subsequent lift-off process. 
Finally, the bridge is underetched in an isotropic oxygen 
plasma which reduces the height of the polyimide layer 
by about 300 to 500 nm. At the narrow constriction the 
nanobridge is now unsupported. The resistance of the 
nanobridge does not increase during the etching process, 
indicating that reaction between the oxygen ions and Zn 
is weak. 

The bridge is then mounted on a three point bending 
mechanism on a stick which is immersed into liquid he- 
lium and can be pumped down to reach temperatures 
of 1.5 K. The residual resistivity ratio of the bridges is 
typically RRR = J R(300K)/i?(1.5K) ~ 2.4 from which 
we deduce an elastic mean free path of 13 nm. On an 
unpatterned thin film of the same thickness we mea- 
sured the critical temperature T c — 0.78 K and a per- 
pendicular critical field of B c = 12 mT. T c is slightly 
reduced with respect to the bulk critical temperature of 
T c b u ik — 0.89 K while B c is enhanced, as usual in thin 
films (Sctaik = 5 mT) (Refill). 

As a first step we have determined the preferred con- 
ductance values of few-atom contacts when opening or 
closing continuously and have constructed the corre- 
sponding conductance histograms. 

A screw with 500 /im pitch, driven by a dc-motor on 
top of the cryostat through an (exchangeable) reduction 
gear box, controls the motion of the pushing rod that 
bends the substrate. The relative displacement of the 
pushing rod can be controlled to a precision of « 4 /nm, 
which due to the geometry of the bending mechanism 
results in a relative motion of the two anchor points of 
the bridge of around 0.9 pm. This was verified using the 
exponential dependence of the conductance on the inter- 
electrode distance in the tunnel regime. The absolute 
precision is restricted due to hysteresis of the bending 
mechanism (torsion of the pushing rod etc.) and rough- 
ness of the screw pitch to about 100 nm for the push- 
ing rod (corresponding to 0.02 nm for the bridge). We 
note that the stability of the setup allows us to halt at 
any point along the breaking or closing curves, even on 
the last plateau when the contact is formed by a sin- 
gle atom, without any observable drift of the measured 
conductance. 

The resistance of the sample is recorded by a 4-point 
resistance bridge with a measuring current of 1 nA and a 
resolution of 10~ 3 in a range from 200 Vl to 200 kfi. Due 
to the mentioned texture of the evaporated film there are 
only few sliding planes parallel to the substrate which 
would allow for a stretching of the film. Consequently it 
often happens that the two electrodes remain separated 
by a gap of several nm after the first breaking. There- 
fore the first opening of the contact has to be performed 
very slowly. Typical opening speeds for the systematic 
measurements are 155 nm/s for the pushing rod, corre- 
sponding to a displacement of the bridge anchor points of 
36 pm/s. The output of the resistance bridge is recorded 
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FIG. 2: a (b and c): Conductance histograms for opening 
and closing the contacts, recorded on sample No. 1 (4) at 
T = 1.5 K. The histograms have been calculated from ~ 
600 individual opening and closing curves, respectively, d (e) 
Examples for typical closing (opening) curves. 



by an oscilloscope and then transferred to a computer for 
calculating the histogram. 

The conductance G decreases in steps of the order of 
0.2 to lGn, their exact sequence changing from opening 
to opening (sec panel c of Fig. 0). Since the mechani- 
cal deformation is different for opening and closing, we 
construct separate histograms for both experiments. Fig. 
12^. displays the conductance histograms calculated from 
~ 600 subsequent opening and closing sequences recorded 
during two weeks on sample No. 1 which has been kept 
at low temperatures. Similar histograms were achieved 
with a lower motion speed of the bridge arms of 17 pm/s. 
After an accidental warming up to w 80 K no step curves 
could be observed any more. 

Several peaks appear for both motion directions at 
similar positions of 0.8 Go, 1.4 Go and w 2.5 Go- In- 
terestingly the first peak is split into two substructures 
at w 0.7 Go and 0.9 Go, a feature which has not been 
observed for other materials. This splitting has been ob- 
served in several independent measurements carried out 
with two different samples (Nos. 1 and 4). Although 
the splitting is observed in all our measurements, the ab- 
solute values of the peak positions and the relative peak 
heights vary slightly from experiment to experiment (see 
Fig. and c) . A possible explanation of the splitting 
would be two different configurations of the central atom 
e.g. different crystalline orientations or a short contact 
with only one plane with a single atom in cross-section 
vs. a dimer contact with two single atoms in series. The 
latter configuration has been found in molecular dynam- 
ics calculations to be typical for Au contacts just before 
breaking 15 . This geometrical difference gives rise to dif- 
ferent conductances^. In similar experiments by Yanson 
et aliii on "classical" breakj unctions made of notched Zn 
wires, no splitting of the first peak is found, but a shoul- 
der that is compatible with our findings. The peaks at 



higher conductance values are in reasonable agreement 
with our findings as well. It is well known that details 
of the histograms do depend on the experimental condi- 
tions, as one can easily see by comparing different pub- 
lished results of e.g. the most extensively studied metal 
gold. In the experiment with Zn by Yanson et al. higher 
voltages have been applied, a slightly higher tempera- 
ture and faster opening speeds have been used and only 
opening traces were recorded. 

In our experiment, the main difference between open- 
ing and closing consists in the peak heights. The low- 
conductance peaks are less well pronounced when clos- 
ing. In the closing traces (see panel d) we often observe 
that the first contact has a conductance well above Go- 
We attribute this to the formation of contacts larger than 
a single-atom contact since the apex atoms from the pre- 
vious opening sequence might have relaxed back to posi- 
tions deeper in the tip. 

Conductance histograms of multivalent metals in gen- 
eral do not show well pronounced multiple peaksiiiS*^. 
An exception of this general observation was so far Al, 
as mentioned already in the introduction. The histogram 
of Zn shows similarities to that of Al: Several peaks oc- 
cur with an average spacing between the peaks of about 
1.3 Go- One could argue whether this similarity in the 
preferred conductance values would be mirrored in the 
quantum transport properties, i.e. the conduction chan- 
nels. Therefore we performed measurements in the su- 
perconducting state for the determination of the chan- 
nel ensemble. For this purpose extremely stable contacts 
and thus a more complex bending mechanism is required. 
The samples are mounted onto a three-point bending 
mechanism thermally anchored to the mixing chamber 
of a dilution refrigerator and that contains now a dif- 
ferential screw close to the sample holder. The pushing 
rod of the bending mechanism is driven by a dc motor 
at room temperature via a mechanical feedthrough and 
several gear boxes with a reduction ratio of 1:2,400,000 
between the turning speed of the motor and the axis en- 
tering the differential screw. The differential screw trans- 
forms one turn of the axis into a motion of the pushing 
rod of 100 fim. The sample geometry is the same as in the 
previous measurements, giving again rice to a reduction 
ratio of 1:5000. Taking into account the reduction of the 
gear boxes, the differential screw and the geometry of the 
sample itself, one turn of the motor yields a displacement 
of the bridge arms of the sample of about 8 fm. There- 
fore, stable tunnel or few-atom wide contacts are achiev- 
able, but no histograms can be recorded within a rea- 
sonable time. For further experimental details see Refs£ 
andSi. Contrary to the experiments described in 5 -, only 
lossy coaxial wires and twisted-pair wires in capillaries 
were used for high-frequency filtering. This restricts the 
electronic temperature of the samples to about 160 mK. 

Fig. depicts a current- voltage sweep (IV) measured 
at T = 160 mK when sample No. 2 was in the close tun- 
nel regime with a conductance of G = 0.34 Go- In the 
inset, the differential conductance is shown. The tun- 
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FIG. 3: Current-voltage characteristic (IV) of sample No. 
2 when in the near tunneling regime at T = 140 mK. From 
the slope at voltages higher than the superconducting gap we 
deduce a conductance of 0.34 Go. The superconducting gap 
is A ~ 165 jj,eV. Inset: Numeric derivative of the IV. 



nel contact was achieved by carefully closing the contact 
right before the so-called "jump to contact", i.e. the 
strong conductance increase signalling htc formation of 
an atomic bond between the two bridge arms 2 ". Appar- 
ently, the spectrum deviates from the expected behavior 
for BCS superconductors, as has already been observed 
in earlier studies of tunnel contacts of Zn fabricated by 
different methods^. We deduce a superconducting gap 
A = 165 jUeV, somewhat higher than the values found in 
the literature, which vary typically between A = 120 /icV 
and 145 /xeV, depending on the measuring method, crys- 
tallographic orientation and other parameters^. After 
these characterization measurements, we close the bridge 
again to form single- or few-atom contacts. This is done 
by reducing the bending of the substrate until the sam- 
ple achieves a conductance of w 10 Go and opening it 
again carefully. Just before break typically conductances 
of 0.7 — 1.4 Go are observed. 

In Fig. 01 we depict a typical set of IVs recorded on 
contacts on last plateaus of sample No. 2 before break- 
ing of the wire, measured at temperatures around 160 
mK. The typical nonlinearities due to multiple Andreev 
reflection (MAR) are observable 22,23 . However, they ap- 
pear rounded compared to the theoretical prediction and 
the observations for contacts of superconductors with a 
density of states obtained from dl/dV vs. V that is 
closer to the BCS shape. This fact limits the precision 
with which the channel ensemble can be determined. The 
lines shown in Fig. 0] are fits to the theory of MAR tak- 
ing into account 1, 2 or 4 channels for curves a, b, and c 
(all measured on sample No. 2). Because of the rounded 
MAR structures, only the number of channels can be es- 
timated, but no meaningful value for the transmissions 
can be given for contacts with more than one channel 
because several channel transmissions combinations give 
fits with the same quality. The fact that the determi- 
nation of the number of channels is more robust than 
the precision with which the individual transmissions can 
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FIG. 4: Current-voltage characteristics of few-atom con- 
tacts of sample No. 2. The measuring temperature was 
T ~ 160 mK. The symbols are experimental data, while the 
fits are calculated with the theory of MAR from Refs. |22Tl23| 
for A = 165 /neV and T = 160 mK. The transmission coeffi- 
cients are: Curve (a) n = 0.97, r 2 = 0.03; (b) n = 0.84; (c) 
ti = 0.69, t 2 = 0.39, r 3 = 0.23, r 4 = 0.15. 



be deduced has already been reported for few-atom con- 
tacts of Nb, a metal that also shows a DOS deviating 
from the BCS shaped. The robustness of the determina- 
tion of the number of channels is due to the super- linear 
dependence of the excess current (the extrapolated in- 
tersection of the high-bias current with the current axis) 
on the transmission values^^ 3 -. Interestingly, contacts 
with a dominating single channel, accompanied by a sec- 
ond, but much smaller one, similar to the observations 
for Au£i24, are achievable. These IVs are routinely ob- 
served at the end of a last plateau after an abrupt change 
of the IV including an increase of the conductance, and 
thus an re-arrangement of the atomic-size contact. Once 
such a contact has been established, it remains stable 
for an elongation of the bridge of 1-2 A. For Al samples, 
contacts with a dominating channel with the contribu- 
tion of one or two smaller channels could be observed as 
wel&iS. However, for Al these contacts evolved continu- 
ously from a situation with three channels of similar size 
and occurred only very rarely^. 

The origin of these single-channel like contacts and 
their possible correlation with the observed splitting of 
the first maximum in the conductance histogram has 
to be clarified further. Therefore, measurements on 
contacts, fabricated from whisker breakj unctions have 
been carried out. The results of these investigations 
together with a detailed theoretical analysis will be re- 
ported elsewherei&. We expect the superconducting DOS 
to be closer to the BCS DOS for those samples and con- 
sequently allowing to determine the transmission values 
precisely. 

Another observation hampers the quantitative deter- 
mination of {t„}. As can be seen in curve c of Fig. 
01 the current onsets corresponding to MAR processes 
involving two or more electrons seem to be shifted to 
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higher voltages with respect to the current onsets pre- 
dicted by the theory and observed for contacts of other 
superconductors^. We attribute this observation to the 
onset of Coulomb blockade^ due to the granularity of 
the film. Our findings are in qualitative agreement with 
a prediction by Avishai et al. 1 . The authors have cal- 
culated the IV characteristics of an interacting quantum 
dot between two superconductors and found a decrease of 
the subgap current and a shift of the higher-order MAR 
processes to higher voltages. However, since no quantita- 
tive theory for our experimental realization and resulting 
transmission coefficients are available, we cannot perform 
a quantitative comparison with this suggestion. 

Cron et aliS 7 . have placed an Al breakjunction in a 
high-impedance environment in order to study the dy- 
namic Coulomb blockade^ of a single junction. They 
found a reduction of the low-bias conductance in reason- 
able agreement with a theoretical model when the con- 
tact was driven normalconducting by applying a small 
magnetic field. However, no effect was detected in the 
superconducting state. 

In order to test the interpretation of Coulomb block- 
ade, we record IVs of tunnel contacts for various mag- 
netic fields, applied perpendicular to the film plane and 
for temperatures below and above T c . Fig. [S] displays a 
typical set of IVs, recorded on sample No. 3 for a contact 
with a conductance at high bias voltages of about 0.3 Go- 
At the lowest measuring temperature and without mag- 
netic field the conductance is reduced for small voltages 
\V\ < 100 /iV. For larger voltages a strong increase is ob- 
served. Around \V\ ~ 200 /xV the conductance decreases 
again. At voltages above 500 /zV the IV is roughly linear. 
The small gap in the IV at low voltage is suppressed by 
applying a magnetic field of only 2 mT, which is smaller 
than the critical field of the unpatterncd film but could 
correspond to the critical field of the lithographiclly pat- 
terned film. We therefore attribute this gap to super- 
conductivity. The remaining nonlinearities would than 
have a different origin. The conductance at intermedi- 
ate voltages is progressively suppressed by increasing the 
field strength. However, the nonlinearities remain up to 
fields of B ~ 4 T, gradually evolving from an enhanced 
conductance at zero bias to a slightly reduced one. 

A possible explanation of the occurrence of such non- 
linear IVs would be a strong enhancement of the critical 
field of a contact with high transmission values r > 0.9. 
A drastically enhanced critical field was observed previ- 
ously for atomic-size contacts of Pb prepared with the 
help of an STM2& in a magnetic field applied parallel to 
the current direction. It was attributed to the incom- 
plete penetration of the magnetic field in narrow struc- 
tures. In our previous measurements on Al contacts and 
Au contacts with Al leads fabricated with the same tech- 
nique and field direction as the present experiments on 
Zn contacts, all signatures of superconductivity vanished 
at magnetic field values smaller than or equal to the bulk 
critical field of Al. However, in those samples the film 
structure was optically smooth. Accordingly, the granu- 
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FIG. 5: Current-voltage characteristics of a tunnel contact 
of sample No. 3 with a conductance of ~ 0.3 Go for the 
magnetic fields and temperatures indicated in the legend. The 
magnetic field was applied perpendicular to the film plane. 
For the curves at B = and B = 2 T data for increasing 
and decreasing voltage are shown in order to demonstrate the 
reproducibility. The nonlinearities observed at low T and low 
B vanish only at T > 1.4 K. 
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FIG. 6: Differential conductance, calculated numerically 
from the IVs of Fig. 5. For clarity, the reproducibility is only 
shown for the curves recorded at B = and B = 2 T. 



larity of the Zn film could be the reason for the unusually 
high critical, the nonlinearities only vanish at tempera- 
tures T > 1.4 K, well above the critical temperature, 
which we estimate to be about 750 mK, since the signa- 
tures of MAR for atomic-contacts realized with the same 
sample disappear at this temperature. 

Therefore and in accordance with the observed shift 
of the MAR peaks we attribute the nonlinearities to 
electron-interaction effects, such as Coulomb blockade, 
resonant tunneling or quantum interference effects. The 
interplay between electron-electron interaction effects 
and superconductivity in atomic-size structures has so far 
not been explored theoretically. Consequently, a detailed 
interpretation of the observed IVs below T c including the 
magnetic-field behavior is not possible. 

Taking the typical voltage scale of the nonlinearities 
of |V| ~ 200 /iV as a measure for the charging energy 
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eV = e 2 /2C, we deduce a capacity C = 0.4 fF corre- 
sponding to the capacity of a sphere with a diameter of 
7 juni. This value is about a factor 10 larger than the 
typical grain size of the Zn film. Thus, we interpret our 
observations to indicate the onset of the interaction ef- 
fect. Due to the low-impedance environment of the cir- 
cuit and the effective screening of the surrounding metal 
all interaction effects would be expected to be damped. 

The numerical derivatives dl/dV of the IVs of FigO 
depicted in Fig. Elshow small, but reproducible fluctua- 
tions superimposed over the larger scale variations. The 
fluctuations have a typical voltage period of V c ~ 30 fiV, 
appear at the same voltage values for both sweep di- 
rections (but different for opposite sign of the voltage) 
and are correlated from curve to curve. These findings 
underline at first that the fluctuations are not due to 
experimental noise and second that the contact has not 
been changed by varying temperature and magnetic field. 
There are several possible origins of the fluctuations, 
including (universal) conductance fluctuations (CF) or 
again the onset of Coulomb blockade due to the granular 
structure of the film. The analysis in terms of CF yields 
a typical coherence length of L c = ^JhD/eV c ~ 800 nm 
roughly corresponding to the grain size of the evapo- 
rated film. A gradual transition from CF to Coulomb- 



blockade behavior has recently been observed in carbon 
nanotubes22*2£. However, the physical mechanism giving 
rise to this transition is not yet clear. 

In conclusion, a single-atom contact of Zn is likely to 
have a conductance of « 0.8 Go- The histograms for both 
closing and opening the contacts indicate the importance 
of different contact geometries. In the IV characteristics 
in the superconducting state we find indications of inter- 
action effects, perhaps due to the granular structure of 
the evaporated films. However, due to the lack of sys- 
tematic measurements and theoretical investigations the 
interpretation of the features being due to a Coulomb- 
blockade type mechanism remains somewhat speculative. 
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